3-junction SQUID rocking ratchet 
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We investigate 3-junction SQUIDs which show voltage rectification if biased with an ac current 
drive with zero mean value. The Josephson phase across the SQUID experiences an effective ratchet 
potential, and the device acts as an efficient rocking ratchet, as demonstrated experimentally for 
adiabatic and nonadiabatic drive frequencies. For high-frequency drives the rectified voltage is 
quantized due to synchronization of the phase dynamics with the external drive. The experimental 
data are in excellent agreement with numerical simulations including thermal fluctuations. 

PACS numbers: 05.40.-a, 05.60.-k, 74.40.+k, 85.25.Dq 



During the last decade, directed molecular motion in 
the absence of a directed net driving force or tempera- 
ture gradient in biological systems has drawn much at- 
tention to Brownian motors. Q, Nonequilibrium fluc- 
tuations can induce e.g. transport of particles along pe- 
riodic structures which lack reflection symmetry - so- 
called ratchets. An important class of ratchets is given by 
the rocking ratchet, characterized by a time-independent 
potential and an external perturbation (driving force) 
which may be either deterministic or stochastic, or a 
combination of both.[E 0, IE IE The number of 
ratchet systems considered for experimental studies has 
steadily been growing during recent years. In par- 
ticular, superconducting ratchets, based on the mo- 
tion of Abrikosov vortic es H , |E El Ill[l3, Josephson 
vortices 0, Q, 0, 0, Il7| or the phase difference of 
the superconduc ting wavefunction (Josephson phase) in 
SQUID ratchets[ll [H HJ H3 have been investigated. 
Those systems offer the advantage of (i) good experimen- 
tal control over externally applied driving forces (here: 
currents), (ii) easy detection of directed motion, which 
creates a dc voltage, and (iii) experimental access to 
studies over a wide frequency range of external pertur- 
bations (adiabatic and non-adiabatic regime) , and transi- 
tion from overdamped to underdamped dynamic regimes, 
enabling studies of inertial effects and transition to chaos. 

Zapata et al. 18] proposed a 3-junction (3JJ) SQUID 
ratchet, which consists of a superconducting loop of in- 
ductance L, intersected by one Josephson junction in one 
arm and by two Josephson junctions connected in series 
in the other arm. For vanishing L a quasi-one dimen- 
sional (ID) ratchet potential can be obtained, and recti- 
fication of an ac bias current for low- and high-frequency 
drive of such a rocking ratchet has been predicted. In 
this paper we investigate a 3JJ SQUID, similar to the 
one proposed in|ll|. We derive the equations of mo- 
tion and the ratchet potential for our type of device, we 
present its experimental realization, and we investigate 
its operation as a rocking ratchet for both adiabatic and 
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FIG. 1: Three-junction SQUID: Equivalent circuit (a) and 
schematic layout (b). 



non-adiabatic drive. We also compare experimental re- 
sults with numerical simulations for our device and for 
the originally proposed 3JJ SQUID ratchet. 

We first discuss the underlying dynamic equations, 
using the resistively and capacitively shunted junction 
modeljlE H3- To simplify the layout, junctions 1,2 in 
the left arm are shunted by a common resistor R s , in con- 
trast to the original proposal [l8j. where junctions 1,2 are 
shunted individually. With Kirchhoff 's laws, the Joseph- 
son equations and the phase differences Sk across junction 
k = 1,2, the current /; through the left arm is 

- + J=It = ^— - S k + ^° (S\ + S 2 ) + Iok sin S k + Ins- 

2 Z7T ilXtis 

( la ) 

$0 = h/2e is the flux quantum; Jjv s is the Nyquist noise 
current (with spectral density Si(f) — AksT/Rg) pro- 
duced by the shunt. We neglect the (subgap) resistances 
of the unshunted junctions which are much larger than 
R s . Note that the original model^] can be obtained 
from ((Taj) by replacing $ (<5i + <5 2 )/27r.R s by $o<4/27i\R fc . 
The current I r through the right arm is 

-- J = I r = ^f^-S 3 + -^-4 + ^03Sin(53 + /Ar3 . (lb) 
I Z7T 27TiX3 

The total bias current is J = Ii + I r f and the Langevin 
equations fTajl . (|lb|) are coupled via the circulating cur- 
rent around the loop J = (/; — I r )/2. The phase differ- 
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ences 6k are connected via 

6 3 -(6i+6 2 ) = ^($ a + LJ) = ^ T . (2) 

The total flux <3>t through the SQUID loop has contribu- 
tions from the applied flux $ a and from the circulating 
current J. For comparison with experimental results, 
we performed numerical simulations to solve the coupled 
Langevin equations (jlap. (|lbl) and J3J). From the sec- 
ond Josephson relation we obtain the momentary volt- 
age across the junctions (Uk = $o^fc/27r) and, by time 
averaging, the dc voltage V across the SQUID and thus 
the current voltage characteristic (IVC) and the critical 
current I c vs. applied flux $ a . 

In the following we derive the ratchet potential in the 
overdamped limit, i.e. with 0ck = 27r/ofc-R|C/b *C 1 for 
all three junctions. The displacement currents oc C^Sk 
in Ijlafl . (jlbjl can then be neglected and the equations of 
motion reduce to only two differential equations 



Ir = 



$0 
2ttR : 

_$o_ 

2irR : 



(Si + 62) +ioisin<5i + I Ns , (3a) 
fa + I03 sin $3 + !n3 , (3b) 



plus conservation of the supercurrent in the left arm 
I01 sin 61 = I02 sin^ 2 , which can be used to express the 
overall phase difference for the left arm Si = Si + S2 in 
terms of Si only 



Si = Si + arcsin (q sin 5\ ) , 



(4) 



with q — Iqi / Iq2 . Without loss of generality we assume 
I01 < I02, i-e. q < 1. Equation can be reversed to 



Si = arctan f — — — - — \ 
Vcoso; + q) 



(5) 



Inserting J^} and the expression for the circulating cur- 
rent J2J into l|3bjl and l|3a|) leads to 



2^R 3 63 - 2 +lN3 



-I 03 sin<5 3 

dU 
~dSl> 



_ { S 3 -Si-2n-) 



$0 j I T T . , . ( sin Si 

„ Oi -- + Ins = -hi sm (arctan — — 

2nR s 2 Vcosd; + q 



(6a) 
(6b) 



+ 2^ {S3 - Sl - 27r *- } = -d5i- 
The potential U(Si, S3) is obtained by integration as 

U(S h S 3 ) = -sloi cos 6 3 + P-(6 3 -S, -2it% 2 

Ty(v / l + 2gcos5 / + g 2 +g-l) . (7) 

The minus (plus) sign in the second line applies for 
cos(<5;/2) > (< 0). The symmetry of U is determined 




FIG. 2: Potential U from eq.© for s — 1/2, and two values 
of applied flux $ a . 



by the applied flux and s = Io 3 /Iqi- If the screening pa- 
rameter (3 L = 2L/ /$ < 1 (with I = (io3 + ^oi)/2), 
the phase differences are fixed by the applied flux $ a ie. 
^3 ~ Si — 27r$ Q /$o; hence ^3 can be eliminated: 

U(Si) = -sI 01 cos(Si + 2tt^) 

+ — (\/l + 2qcosSi + q 2 + q - 1) . (8) 

If the two junctions in the left arm have equal critical 
currents (q = 1) the above expression turns into 

U(Si) = -s/01 cos (Si + 2^) - 2/oi cos (Si/2) , (9) 

which coincides with the result presented in^l and is 
shown in Fig. [3 for s = 1/2 (the value proposed infill]). 

As one can see from Fig. [21 it is the applied flux that 
breaks the reflection symmetry of the potential. The crit- 
ical currents /+(<E> a ), J~($ a ) represent the force to over- 
come the potential barrier for positive and negative bias 
current polarity, respectively. Hence, measuring if- pro- 
vides a way to probe the asymmetry of the potential. The 
maximum of I c ($> a ) — 2io is reached if each arm carries a 
current, that is equal to its critical current, i.e. I r = Iq 3 
and = Iqi. The circulating current is then given by 
J = (Iqi ~ -^03)/2, and the three phase differences have 
the values Si — S 3 = n/2 and 62 = arcsin q < ir/2. In- 
serting this into J3J) and solving for $ gives the applied 
flux at maximum I c (+ and — corresponds to positive 
and negative bias current, respectively) 

^ = (- arcsin q + 2(3 L i^i) . (10) 

4 7T 1 + 8 



If q w 1, the first term yields a shift f» =F$o/4 of the 
7^ = ($ a )-dependence. The second term has the same sign 
and is proportional to (3l- Equation ljTU|) can be used to 
estimate q, 0l and s from the measured I c ($ a ). 

The 3JJ SQUIDs we investigated were fabricated at 
HYP RES [2j| using Nb/Al-AlO^/Nb junctions of nomi- 
nal critical current density jo — lkA/cm 2 at T=4.2K, 
capacitance per area C — 38fF//im 2 and nominally 
identical shunt resistors R s — R 3 . Junction areas Aj i — 
Aj t 2 — 3.8 x 3.8^m 2 , and Aj t3 — 2.6 x 2.6/im 2 , corre- 
spond to s w 1/2 and q = 1. For microwave irradiation 
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(up to 28 GHz), the SQUIDs are integrated in a copla- 
nar waveguide (see Fig. ^b)) with 50 SI impedance. In 
total we investigated six devices which differed by the 
size of the SQUID hole, i.e. by the SQUID inductance 
L. All devices showed very similar behavior. Below we 
discuss only one device with a 2 x 3/im 2 hole. All mea- 
surements were performed at T = 4.2 K in a magnetically 
and electrically shielded environment, with low pass fil- 
ters in the voltage leads. Apart from microwave drives 
we also studied the ratchet behavior in the kHz regime 
with ac currents directly fed into the bias leads. 

Fig. [St a) shows a measured IVC (<£ a = 0) with normal 
resistance R n w 0.7 SI, and thus R s w R3 = R « 1.4 fl. 
The inset shows the critical currents J^($ a ) with max- 
imum I c = 2I = 186/iA, giving a noise parameter 
T = 2irk B T I (Jo$o) ~ 2 • 10~ 3 , characteristic voltage 
V c = IqR — 130 ^iV and characteristic frequency f c = 
Vc/^q = 63 GHz. With the design values for the junc- 
tion capacitances C\ = C2 = 0.54 pF and C3 = 0.27 pF 
we obtain /3 c ,i = Pc,2 = 0.38 and /3 C ,3 = 0.09. The 
shape (asymmetry) and modulation depth of / c ($ a ) is 
determined by s, q and (3l- Adjusting these parame- 
ters in numerical simulations of / c ($ a ) (with values for 
r and Pck as determined above) , we find best agreement 
with the measured 7 c ($ a ) for (3^ = 0.1, s = 0.5 and 
q = 0.99. These parameters are very close to the design 
values, and have been used later to numerically calculate 
the response of the SQUID to an ac bias current. 

As a prominent feature in the experimental zero field 
IVCs, we observe step-like structures in the resistive 
state; simulated IVCs exhibit the same generic shape 
[cf. Fig. Of a)]. The simulation allows to trace the in- 
dividual voltage drops across the junctions, as shown in 
Fig. injb). It turns out that the step structures in the 
IVCs can be associated with switching of junctions 1 and 
2 between different states: (i) one junction is in the zero 
voltage state while the other carries a large voltage; (ii) 
both junctions are in the voltage state and oscillate co- 
herently, i.e. with the same frequency; typically these 
oscillations are out of phase [cf. upper inset in Fig.^b)]; 
(iii) both junctions are in the voltage state; however, os- 
cillations are irregular and incoherent [cf. lower inset in 
Fig. Efb)]. This state is associated with a more "noisy" 
IVC as compared to (ii). Our simulations show that over 
a wide range of bias currents the junctions oscillate in- 
coherently and thus are far from the scenario discussed 
in jl8l| . where junctions 1 and 2 were assumed to behave 
identical, i.e. Si = 82 at all times ("synchronous state"). 

We performed extensive numerical simulations to com- 
pare our device (with common shunt for junctions 1,2) 
with the one in |li|. For T = both devices exhibit 
phase synchronous oscillations over the entire range of 
bias currents if junctions 1,2 are identical and with ini- 
tial values 5k — and 8k — 0. However, choosing either 
different initial conditions or non-identical junction pa- 
rameters results in a behavior very similar to the one 
shown in Fig. Ofb). At T = the synchronous state 
is more easily destroyed if junctions 1,2 have a common 
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FIG. 3: SQUID characteristics without ac drive: (a) I(V)- 
curves (<3? a = 0) and I c (&a) (inset), (b) Simulated V(J)-curve 
(black) and individual voltages across junction 1 (grey) and 2 
(light grey). Insets: Voltage vs. time for junction 3 (black), 1 
(grey) and 2 (light grey) for I/Iq = 3.2 and 2.8. Parameters 
for numerical simulations: T = 2 • 10 -3 , /3c,i = Pc,2 = 0.38, 
Pc,3 = 0.09, Pl = 0.1, s = 0.5 and q = 0.99.' 



shunt, as compared to being shunted individually. In the 
latter case, the synchronous state is more stable, due to 
a damping term cx d(5i — 82) /dt which is absent in the 
equations of motion for the system with common shunt 
for junctions 1,2|24|. However, thermal fluctuations at 
finite temperature, also tend to destroy the synchronous 
state for the individually shunted 3JJ SQUID. Hence, for 
the experimentally relevant case, i.e. at finite tempera- 
ture (r^lO -3 ) and slight junction asymmetries (q ^ 0), 
both types of SQUID ratchets show typically incoherent 
oscillations of junctions 1,2. 

Fig. 0] shows experimental data, and numerical simu- 
lation results, for the rectified voltage V under ac drive 
/ = I ac sin 2wft, clearly demonstrating operation of the 
SQUID as a rocking ratchet. For adiabatic drive / <C f c 
[cf. Fig.QJa)] rectification appears at I ac ~Io- V sharply 
peaks at I ac w 2Iq and decreases again with increasing 
I ac , as predicted for a ID rocking ratchet |j| and as ob- 
served experimentally for asymmetric (two junction) dc 
SQUID ratchetsGjlllJ. In contrast to the prediction for 
overdamped rocking ratchets 0, Q , the initial increase in 
V(Iac) up to the maximum V max is not linear. Instead, 
it sets off with a very steep slope dV/ dl ac which becomes 
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FIG. 4: Rectified voltage i«. ac drive amplitude at <& a — 
$o/4 for a 3JJ SQUID ratchet with increasing drive frequency 
/ (from top to bottom). Left column shows experimental 
data at T = 4.2 K (insets in (c),(e) show enlarged section of 
one voltage step) and right column shows the corresponding 
simulated curves (parameters as in Fig.^J. The dashed lines 
in (b), (d), (f) were calculated for T = 0, with Joi = Joa 
(q = 1) and 8 k = 5 k = (fe = 1, 2) at t = . 



smaller at higher voltage. A strikingly similar behavior 
has beenpredicted for a strongly underdamped rocking 
ratchet 7J. In our case, the motion of Si is overdamped; 
however for 5i — 62 it is underdamped, which obviously 
gives the same result as reported in [3. Vm ax = 25/xV is 
as large as 20% of V c , in excellent agreement with sim- 
ulations shown in Fig. Efb) . As maximum rectification 
requires $+ — $~ w $ /2, it is easy to understand from 
ca. (|10|l why the 3JJ SQUID ratchet is superior to the 
asymmetric dc SQUID ratchet: For the asymmetric dc 
SQUID the arcsin-term in (|10|l is absent, and the shift 
in the J^($ a )-curves is determined by [3l and s only. 
Hence, the condition for maximum rectification depends 
strongly on the SQUID parameters and is more difficult 
to fulfill experimentally. In contrast, for the 3JJ SQUID, 
it is sufficient to keep (3l small, so that the arcsin-term 
dominates. For q « 1, the condition for maximum recti- 
fication is then easy to achieve experimentally. 

At nonadiabatic drive frequencies [c.f. Fig.^c),(e)] we 
observe a step-like increase of V(I ac ) up to V max and 
oscillations in V(I ac ) for higher drive amplitudes. The 



voltage steps appear at V/(IqR) — n ■ (f/f c ) {n: inte- 
ger) and can be interpreted as Shapiro steps, where the 
phase dynamics synchronizes with the external drive pcj. 
This quantization of the ratchet effect has been predicted 
as a characteristic feature of nonadiabatically driven ID 
rocking ratchets 0, however, it has not been ob- 
served experimentally until now, due to thermal noise 
smearing |20j. For the same reason, V(I ac ) oscillates at 
higher drive amplitudes, rather than showing a step-like 
behavior. This effect of thermal noise was also predicted 
in 4] for ID rocking ratchets and is clearly demonstrated 
for our device in the simulations of V(I ac ) shown in 
Fig- Eld), (f) for finite and zero T. For finite T, our sim- 
ulations nicely reproduce the experimental data. For the 
case T = 0, simulations were performed with q — 1 and 
identical initial conditions 5\p — £1,2 = 0. This leads to 
the formation of the synchronous state (i.e. Si — 5% at all 
times). Only in this case a step-like switching of the rec- 
tification is observed in the numerical simulations also at 
large ac drive amplitudes; i.e. experimental observation 
of quantized rectification for large driving amplitudes, 
will be quite difficult. Furthermore, we note that only in 
the synchronous state {(f = 81 — 82 = 0) the dynamics 
of the system is described by Si(t) only, as the motion 
of a particle in the 1-dimensional potential 10 which is 
47r-periodic in <5/. In contrast, if <fi ^ 0, the dynamics 
can be described in a 2-dimensional potential along Si 
and ip [24j. This potential has minima which are sepa- 
rated by 2n if projected onto the 5/-axis. This explains 
the doubling of the Shapiro step height calculated for the 
synchronous state [c.f. dashed lines in Figs. EJd),(f)] as 
compared to the calculations for finite T and as observed 
experimentally [c.f. FigsQJc),(e)]. 

In conclusion, following a proposal in|l8|. we have re- 
alized 3 JJ SQUIDs and demonstrated operation of those 
devices as very efficient rocking ratchets. In contrast to 
the original proposal, which assumed that the two junc- 
tions connected in series oscillate in-phase, we found that 
these junctions usually oscillate incoherently. Nonethe- 
less the devices show a large ratchet effect, more than two 
orders of magnitude above initial estimates given in|l8|. 
Apart from rectification of a harmonic low-frequency 
drive, we find rectification of nonadiabatic drives with 
striking properties, such as Shapiro- like steps, i.e. quan- 
tization of the velocity of directed motion, which depends 
only on the drive frequency. The simplicity of the de- 
sign, the good experimental control over external drives 
and device parameters which define e.g. the shape of 
the ratchet potential, and, finally, easy detection of di- 
rected motion by measuring voltage, gives this device ex- 
cellent perspectives for further basic experimental studies 
of ratchet effects, e.g. for stochastic drives. 
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